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Abstract 

The high-T c cuprate superconductors are now believed to be intrinsically inhomogeneous. We 
develop a theory to describe how this inhomogeneity affects the intermodulation coefficient of 
such a material. We show that the continuum equations describing intermodulation in a super- 
conducting layer with spatially varying properties are formally equivalent to those describing an 
inhomogeneous dielectric with a nonzero cubic nonlinearity. Using this formal analogy, we cal- 
culate the effect of inhomogeneity on the intermodulation coefficient in a high-T c material, using 
several assumptions about the topology of the layer, and some simple analytical approximations to 
treat the nonlinearity. For some topologies, we find that the intermodulation critical supercurrent 
density Jimd is actually enhanced compared to a homogeneous medium, thereby possibly leading 
to more desirable material properties. We discuss this result in light of recent spatial mappings of 
the superconducting energy gap in BSCCO-2212. 
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I. INTRODUCTION 

Ever since the discovery of high-T c superconductors 1], many workers have attempted 
to develop practical applications for them. One potential electronic application is as a 
microstrip resonator. Such a device has been developed by Willemsen et al. |2j, using a 
high-T c cuprate material. Even though these devices do not work at very high current 
densities, they are subject to strong nonlinear effects which mix together microwaves of 
different frequencies. This mixing, known as intermodulation, was studied experimentally 
and theoretically by Willemsen et al. Q; theoretical models to explain their measurements 
were developed by Dahm and Scalapino N, m . 

n n 

In the model of Refs. [3|] and the intermodulation is described in terms of an intermod- 
ulation current density Jimd, which depends on both temperature and the angle between 
the direction of current flow and the crystal axes. An equivalent quantity was also considered 
by Yip and Sauls for a ci-wave superconductor. Because of low-lying quasiparticles, they 
found that Jimd was much smaller than for a corresponding s-wave superconductor and also 
that it depended on the angle between the in-plane current density and the k-vector of the 
gap nodes. Several other workers have also studied this intermodulation and the harmonic 
generation due to nonlineax effects flQS. 

In all this work, the high-T c cuprate superconductor was assumed to be homogeneous - 
that is, the Cu02 layer properties were assumed independent of position within the layer. 
However, recent experimental work has invalidated t his p icture. Specifically, in optimally 
oxygen-doped underdoped, and slightly overdoped jl(| B^S^CaC^Os+s (usually called 
BSCCO-2212), the superconducting energy gap was found to be spatially varying. This result 
was obtained from scanning tunneling microscopy/spectroscopy (STM/S) images of the 
superconducting layers. The differential tunneling conductance spectra of the sample were 
measured and the position-dependent gap A was inferred from measurements of the energy 
difference between two coherence peaks in the spectra above and below the Fermi level. In 
the underdoped sample, the gap was found to map into two distinct regions. One (called the 
a domains) had a gap A < 50 meV; the other (denoted the (3 domains) had A > 50 meV. 
The a-domains were identified as superconducting, because they showed coherence peaks 
in the tunneling spectra. The /^ -reg ions were found to be non-superconducting, and were 
identified as a pseudogap phase |n|, with a large gap. It was concluded that the spatially 
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varying superconducting energy gaps do not arise from impurities, but are instead intrinsic 
properties of the material. Thus, underdoped BSCCO-2212 is an intrinsically granular 
superconductor. 

In this paper, we consider how this inhomogeneity affects the intermodulation in a high- 
T c superconductor, such as BSCCO-2212. As implied above, this intermodulation coefficient 
is actually a nonlinear transport coefficient. In fact, this coefficient describes the current- 
dependence of the superffuid density in the superconductor. This current-dependence is 
particularly strong in the <i-wave high-T c materials, because quasiparticles are excited even 
at very low applied currents. We will show that Jimd is formally analogous to another 
coefficient, well known in the study of nonlinear dielectrics. Using this connection, we will 
demonstrate that Jimd is very sensitive to the detailed geometry of the superconducting 
inhomogeneity. 

We will consider only two-dimensional (2D) systems, and only very low frequencies. This 
regime is appropriate to the cuprate superconductors, where superconductivity is believed 
to occur within Cu02 planes. Our low-frequency approach should be applicable so long as 
the length scale of the inhomogeneity is much smaller than the wavelength of the applied 
microwave fields, a condition easily satisfied at microwave frequencies. 

The remainder of this paper is organized as follows. In Section |Hj we present the formal- 
ism for calculating the enhancement of Jimd i n an inhomogeneous 2D superconductor. In 
Section HU we give several analytical results for the relevant enhancements and for Jimd-, 
based on the analogy to a nonlinear dielectric composite. Section ITVl presents a concluding 
discussion. 

II. FORMALISM 

A. Intermodulation Coefficients from Ginzburg-Landau Theory 

We begin by expressing Jimd in terms of coefficients of the Ginzburg-Landau free energy 
density F. In the absence of a vector potential A, F takes the standard form 



where a, b, and C are appropriate constants, and ip is the complex position-dependent 
Ginzburg-Landau wave function. For the present problem, we will eventually assume that 




(1) 
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all three constants are functions of position. We also write C = h 2 /(2m*), where m* is a 
quantity with dimensions of mass. 

The local supercurrent density takes the usual quantum-mechanical form 

' [fHfi^ + ac], (2) 



2m* 

where e* is the charge of a Cooper pair and 'c.c.' denotes a complex conjugate. We write 
ip = \ exp(i(f)) and initially assume that <p, but not is position-dependent, so that 

J = 4h/f(W0). (3) 

In the limit of very small current density, the wave function is found by minimizing the free 
energy with respect to \ip\ 2 , which gives the standard expression \ip\ 2 = —a/b. This quantity 
is positive if a < 0. 

If there is a finite phase gradient, F takes the form 

F = aH 2 + ^ + ^H 2 |V0| 2 . (4) 

Minimizing with respect to the modulus of the wave function at fixed V0, we find that 

[q + (ft 2 /(2m*))|V</,| 2 ] 

m = 1 • ( 5 ) 



The corresponding current density takes the form 

2 

(6) 



e 

m* 



In the above discussion, we have assumed that \ is position-independent, so that V|^| = 
0. Even if V|"0| 7^ 0, Eq. (jSJ) for J remains valid. However, there is an extra term in the free 
energy density; so Eq. (JHJ), and hence Eq. (jSJ, do not hold exactly. Nevertheless, we shall 
assume that the most important effects of inhomogeneity can be included by writing 

J = ir 1 (x)V0(l-ir 2 (x)|V0| 2 ), (7) 

with appropriate coefficients -Ki(x) and K 2 ('x). 

We now show that Eq. (jUJ) (for a uniform superconductor) contains the intermodulation 
phenomenon of interest. First, we rewrite Eq. (JBJ) for a uniform superconductor as 

3 = K 1 V<P(1-K 2 \V<P\ 2 ), (8) 



where K\ and K2 are related to the original a, b, and m*. If A 7^ 0, in order for the gradient 
to remain gauge- invariant, we must make the replacement —ihV — > —ihV — (e*/c)A, or 
equivalently 



V0 -> V0 - — A. 

n.c 

Thus, if A 7^ but the phase is uniform, we must rewrite Eq. (JBJ) as 



he 



A 



1 - 



ft 2 c 2 



IAI 



(9) 



(10) 



To order A 3 (or J 3 ), we can replace |A| 2 on the right-hand side of this expression by 
[hc/(e*Ki)] 2 \3\ 2 , which gives 



he 



A 



J 2 



T 2 

J IMD 



where 



K 2 

J IMD — is ■ 



(11) 



(12) 



Finally, we show that Eq. (fTTj) implies a current-dependent penetration depth. To see 
this, we first take the curl of Eq. (jUj) in the low current limit to get 



V x J = -a'B, 



(13) 



where a' = e*Ki/(hc). When this equation is combined with Ampere's law, V x B = 4irJ/c, 
we get V 2 J - [1/A 2 (T, 0)] J = 0, where 

1 47^*^1 



A 2 (T,0)~~ he 2 ' 

and A(T, 0) is the zero-current penetration depth at temperature T. Thus finally 

J 2 1 



(14) 



A 



1 - 



T 2 

J IMD 



(15) 



4ttA 2 (T,0)' 

Eq. (HHD has the form J = -//(T, J) A, where fi(T, J) = c/[4vrA 2 (T, J)], A(T, J) being the 
temperature- and current- dependent penetration depth. Thus, Eq. (fTK|) is equivalent to 



J 2 



T 2 

J IMD 



A 2 (T,J) A 2 (T,0) 
To order J 2 , this result is equivalent to Eq. (7) of Ref. 



(16) 
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B. Estimate of Ginzburg-Landau Parameters for Cuprate Superconductors 



Within the Ginzburg-Landau formalism, the penetration depth is related to the order 
parameter by 

"47re* 2 |^| 2n 1/2 



1 



A(V.O) 1 ' (17! 

and |^| 2 = —a/b, provided that a < 0. Taking a = — a(l — T/T c q), where a > and T c q is 
the mean-field transition temperature, we obtain 



\^\ 2 = a(l-T/T c0 )/b 



and 



1 



47re* 2 a(l - T/T c 



cO ) 



(18) 



(19) 



A 2 (T,0) m*c 2 b 
In Eq. (|19|). -0 which has dimensions of a wave function. Lacking an accepted microscopic 
theory for the high-T c cuprates, we may estimate a and b using BCS theory, as discussed, 
for example, by de Gennes [l2|. The result is 

h 2 



a 



and 



h 4 



0.1066 



4m* 2 ^ N{0){k B T, 



cO 



\2 ■ 



(20) 



(21) 



Here, N(0) is the single-particle density of states at the Fermi energy (measured in states 
per unit energy per unit volume), and £ is a temperature- independent coherence length. 
The penetration depth is then found to be determined by the equation 

1 32ne 2 



wop M-xwiW, 

where A(T) is the equilibrium value of the energy gap, given by the equation 



(22) 



A 2 (T) = 9.38fc|T c0 (T c0 -T). (23) 

Now according to experiments 0, Q], in the small gap regions, BSCCO-2212 has a 
sizable superfluid density, whereas in the large-gap regions, the superfluid density is small 
or zero. If we interpret 1/A 2 (T, 0) as proportional to the local superfluid density, then this 
experimental result implies that 1/A 2 (T, 0) should vary inversely with A(T). In order for 
this to be consistent with Eq. (|52|L the quantity iV(0)£ 2 A 2 (T) should vary inversely with 
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A(T). We therefore assume that 1/A 2 (T, 0) oc A~ X (T), where x > is some number which 
characterizes BSCCO-2212. While this is a highly oversimplified model, it does suggest how 
Jimd is influenced by the inhomogeneity of the high-T c layer. 

III. MODEL CALCULATIONS 

A. Analogy to a Composite Dielectric Medium with a Cubic Nonlinearity 

We now apply the above results to calculate Jimd for several models of inhomogeneous 
superconducting layers. In all cases, we attempt to choose the layer properties to resemble 
those reported in experiments on BSCCO-2212. Our goal is to solve Eq. ((7|) for J(x) and 
0(x) for some prescribed inhomogeneous superconductor. We assume that -ft'i(x) and K 2 (x.) 
are specified, but random. 

The present problem is formally equivalent to a randomly inhomogeneous dielectric with 
a cubic nonlinearity Q, In the latter problem, the electric field E and electric displace- 
ment D are related by 

D(x) = e(x)E(x) + X (x)|E(x)| 2 E(x), 
V x E = 0, 

V-D = 0. (24) 

For the intermodulation problem, the analogous equations are Eq. 0, supplemented by the 
steady-state charge conservation condition V • J = 0. Thus, — V</> plays the role of E in 
the intermodulation problem, and the phase plays the role of the scalar potential. The 
quantities i^i(x) and i^i(x)i^2(x) are analogous to the linear dielectric function e(x) and 
cubic nonlinear susceptibility x( x )- The quantity — is, of course, curl-free like E in an 
electrostatic problem. Thus, we are again connecting a divergence-free field to a curl-free 
field. To treat the intermodulation problem, therefore, we can use all the formal results 
previously derived for an inhomogeneous nonlinear dielectric, which we now briefly review. 

For a material described by Eqs. (j2U), two useful coefficients are the effective linear 
dielectric function e e and the effective cubic nonlinear susceptibility \ e . These quantities are 
defined in terms of the space-averaged electric field (E) and displacement (D) by 

(D)= ee (E) + Xe |(E)| 2 (E). (25) 
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As shown in Ref . ^| , \e can be expressed as an average over the fourth power of the electric 
field in the associated linear composite. That is, if Eii n (x) is the electric field in a composite 
described by the linear relation D(x) = e(x)E(x), then \ e is given by 

Xe £ 4 = ( X (x)|E lm (x)| 4 >, (26) 

where E Q is the applied electric field. If the composite medium has n components, the ith 
of which has nonlinear susceptibility Xi, then Eq. (|2l)j) can be rewritten as 

Xe 



^PiXiCi, (27) 



where Pi is the volume fraction of the ith component, is an enhancement factor given by 

„ (|E(x)| 4 ) iilin 



FA 



(28) 



and (• • ■ }i i ii n means a space-average within the ith component in the related linear medium. 
Thus, Ci describes how much the fourth power of the electric field is enhanced in the ith 
component in the linear limit. 

The moments are difficult to compute exactly, except in a few very simple geometries. 
We have therefore chosen to make a decoupling approximation ^ Q| j specified by 

(|E(x)|% in ^(|E(x)|% in . (29) 

Clearly, the decoupling approximation will be most accurate if the fluctuations 

(|E(x)| 4 ),;.iiTi — (|E(x)| 2 )^ lin within the ith component are small compared to (^(x)! 4 )^ 
itself 16j . This assumption is most likely to be accurate in geometries such that |E(x)| is 
uniform in each of the nonlinear components, but will be less accurate when the fluctuations 
are large. For example, in the so-called Hashin-Shtrikman geometry, in which one of the 
two components is embedded in the other, these fluctuations are small in the embedded 
component, and hence this approximation will be excellent if only the embedded component 
is nonlinear. However, in a composite near a percolation threshold, the fluctuations will be 
large and this approximation will be less accurate. 

If we make the approximation (}29|) . we can express (^(x)! 2 )^ exactly in terms of the 
effective linear dielectric function e e through the relation 

Pi (|E(x)| 2 ) lilin _ de 

J 



Here, <9e e /<9e; is the partial derivative of e e with respect to e», at constant £j (j ^ i) and con- 
stant volume fractions pj. Given a simple analytical approximation for e e , these derivatives 
can be easily computed in closed form, thus yielding a simple analytical approximation for 

Xe as 

Xe = Y^XiF?/Pi (31) 

i 

with ej = Ff /pj. We will use this approach, combined with the analogy described above, to 
obtain approximations for the intermodulation coefficient in an inhomogeneous supercon- 
ducting layer. 

In the present work, we use two different approximation methods to calculate e e : the 
effective-medium approximation (EMA) 0,0], and the Maxwell- Garnett approximation 



(MGA) [18]. The EMA is suitable for a binary composite with symmetrically distributed 
components, so that neither can be viewed as the inclusion or the host 15]. In this case, if 
the components are isotropic, e e satisfies the quadratic equation 

PA " A 7 te ; + (1 - Pa) " B ~ te r = 0. (32) 
e e + g{e A - e e ) e e + g{e B -e e ) 

Here pa is the volume fraction of the component A, €a and e# are the dielectric functions 
of the components A and B, respectively, and g is a 'depolarization factor': g — 1/2 in two 
dimensions (2D) and g = 1/3 in 3D. The physically meaningful solution of Eq. (|32|) is the 



root which varies continuously with pa, approaches the correct 



imits ai pa = and 1, and 



3. 



has a nonnegative imaginary part when and e# are complex 

The MGA is more suitable to a binary composite where one component can be regarded 
as a host in which the other is embedded |l6(. When the host material is isotropic and 



linear, the effective dielectric function of the composite takes the form 

Pi (e< - e h ) 



1 + 



(33) 



;i -pi)[e h (l -g) + eig] + p t t h _ 
where pi is the volume fraction of the inclusion, and are the dielectric functions of the 
inclusion and the host, respectively, and g is again the depolarization factor. 



B. Application to an Inhomogeneous Superconductor 

We can readily use the above analogy to compute the effective nonlinear coefficients for 
an inhomogeneous superconducting layer. We consider a superconducting layer comprised 
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of two "components," A and B, with areal fractions p A and p B = 1 — Pa, which have two 
different energy gaps. The two components are both intrinsic to the given sample, in the 
sense that they are not caused by the introduction of non-superconducting impurities. A 
realistic sample of BSCCO-2212 probably has a continuous distribution of gaps, but we 
make this simplification for computational convenience. 

The effective cubic nonlinear coefficient (K x K 2 ) e takes the form 

(K x K 2 ) e = (K x K 2 ) A p A e A + (K x K 2 ) B p B e B . (34) 

Jimd m Eq. (fT2"j) thus becomes 

K u (K x K 2 ) e 

JlMD ~ ~ ^3/2 • W 

A 2e A 2e 

To apply the present formalism, we need to find suitable K x and K 2 values. From Eqs. 
® and ©, K x = -(a/b)(he* /m*) and K X K 2 = h 3 e* / (2m* 2 b) , giving K 2 = -h 2 /(2m*a). 
Using a = a(t — 1), where t = T/T c0 , and taking a and b from Eqs. (fTU|) and (j20J), we find 

a fie* fie 2 

Kl = ~b~^ = 47re*A 2 (t,0) ^ 36 ^ 

and 

K 2 = = e 2 (t). (37) 

2m*\a\ s w v ; 

In typical cuprate superconductors, £(£ = 0) ~ 15 A and A(t = 0, J = 0) ~ 1500 A; so 

K 2 = (15A) 2 . (38) 

To estimate the values of K x , we first assume that 

1 



oc [A(T)]" X . (39) 



A 2 (T,0) 

This assumption embodies the experimental observation that the superfluid density is rel- 
atively large in regions where the gap is relatively small. Thus, it is simply an effort to 
include relevant experimental features in the model, without attempting to explain them. 
In the model calculations below, we consider two different values of x, in order to see how 
this value affects the calculated Jimd- 

Eqs. (I36J) and (I39|) can be combined with experiment to get a rough estimate of K x . 
According to Ref. |9(, A(T) ranges from 25meV to 65meV in optimally doped BSCCO-2212 
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at low T. We assume that the mean value, 45 meV, corresponds to A(0, 0) = 1500 A. This 
fixes the proportionality constant in Eq. (}3*9"j) . Using this proportionality constant and Eq. 
(I55J) . we get 



K x = 3.49 x 10 11 ( ~^%) (10) 
\Ao meV / cm • s 

We first assume that x = 1. Eq. (|4()jl implies that A(0) = 25meV and 65meV correspond 
respectively to K 1A = 6.28 x 10 11 esu/(cm • s) and K 1B = 2.42 x 10 11 esu/(cm • s). For K 2 a 
and K 2 b, we have no definite information from experiment. We therefore assume simply 
that K 2A = K 2 b- 

In Figs. [T]and|2J we show the calculated enhancement factors and e# and the corre- 
sponding intermodulation critical current density Jimd for these models, as functions of pa- 
In Figs. Ufa), QIc), andEfa) we use the EMA [Eq. @2J], while in Figs. Hfb), OJd), andEfb) 
we use the MGA [Eq. (jSHJ)], with B considered as the host. In both cases, we combine these 
approximations with the decoupling approximation [Eqs. (|29 |) -(|5T |) ]. to obtain Jimd- 

Fig.Efa) shows that Jimd increases linearly with pa in the EMA. As in Eqs. ()34|) and ()35jh 
Jimd has contributions from the nonlinearity of both components. While the enhancement 
factor tA is never larger than unity, e# can exceed unity, depending on the value of pa- 
Thus, the nonlinearity of B has a larger influence on Jimd than that of A. As a result, 
Jimd behaves similarly to eg, having a larger enhancement for the larger Kia/K\b- The 
MGA results differ little from the EMA results except for a broad peak around pa = 0.9 for 
the larger ratio of Kxa/Kxb- This peak results from the shift to higher values of pa of both 
Pb&b and Pa^a, seen in the MGA results of Fig. [TJ 

We can also calculate the effective linear coefficients K\ e for these two models. In the 2D 
EMA, K Xe satisfies 

p A KlA ~ Kle + (1 - PA ) KlB ~ Kle = 0, (41) 
PA K 1A + K le 1 PA) K 1B + K le y ' 

while in the 2D MGA with B considered as the host, we get 

V K [1 -l 2Pa(K 1a -K 1B ) ] 

Kle = I 1 + (1-pa)(K 1 a + K 1b ) + 2paK 1b \ • (42) 
The i^ie's are proportional to the effective superfluid densities (or the effective inverse-square 
penetration depths) of these 2D materials in the linear limit of very small applied currents. 
The values calculated from the EMA and MGA are shown in Figs. Efa) and^b). Both 
increase monotonically with increasing areal fraction of the small-gap component A. The 
MGA results differ very little from the EMA results. 
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In Figs. |H and |HJ we show an analogous set of calculations, but with x = 3. We 
again assume a binary distribution of gaps, using the same gaps as in Figs. UHS1 Because 
of the larger x, the ratio Kia/Kib is larger than in Figs. EtH For x = 3, using the same 
proportionality constant, we find that the gaps A(0) = 65meV and 25meV now correspond 
to K x = 1.16 x 10 11 esu/(cm • s) = K 1B , and K x = 20.35 x 10 11 esu/(cm • s) = K 1A . 

For x = 3, for the larger ratio of Kia/Kib, the enhancement factor pb&b has clear peaks 
as a function of the areal fraction pa- This behavior is shown in Figs. HJc) and^d). The 
peak occurs at around the percolation threshold of pa = 0.5 in the EMA results, but at 
around pa = 0.95 in the MGA results. In addition, pb^b is dramatically larger (~ 100) in 
the MGA than in the EMA, for the larger ratio of K 1A /Ki B . Note also that the EMA results 
are nearly symmetric about pa = 0.5 while the MGA results are very asymmetric. There 
is also a large difference between results for the two gap ratios in the MGA results, but a 
smaller one in the EMA. By contrast, Pa^a is monotonic in either the EMA or the MGA. 
Since Jimd has two contributions, one from the enhancement of A and the other from the 
enhancement of B, one expects that the behavior of Jimd in the MGA results will mirror 
the enhancement factor pb^b because pa&a "C Pb^b for the larger ratio of Kia/Kib- 

Fig. El shows the behavior of Jimd for the K x s shown in Fig. 0] As expected, and as 
already found for x = 1, Jimd for x = 3 generally follows the trend oipB&B- In particular, 
because of the clear peak in Pb^b, the EMA results show a weak broad peak near the 
percolation threshold p c in Jimd for the larger ratio of Kia/Kib- The EMA and MGA 
results differ greatly for the larger ratio of Kia/Kib, not only in the shape of the curves but 
also in the magnitude of Jimd- In this case, the MGA results follow mostly the shape of 
the curve for Pb^b- Although Jimd increases monotonically with pa in the EMA results, it 
drops sharply above pa = 0.95 for the larger ratio of Kia/Kib in the MGA results. Overall, 
the x = 3 case produces a much larger value of Jimd for the larger ratio of Kia/Kib- Thus, 
for a device requiring a large Jimd, these results suggest that the best results would be 
obtained using an inhomogeneous superconductor with a large gap difference and a large x 
in a Maxwell-Garnett geometry. 

Fig. |U1 shows the effective superfluid densities K Xe with x = 3. For the smaller ratio of 
Kia/Kib, the EMA results differ little from the MGA results as we found previously in Fig. 
EJ However, for the larger ratio of Kia/Kib, the two differ significantly. 

There is no distinction between the inclusion and the host in the EMA method, but 
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there is in the MGA method. In our MGA results thus far, we have assumed that A (the 
component with the smaller gap) is the inclusion and that B is the host. We now consider 
the reverse configuration, where B is surrounded by A. Results for this configuration are 
shown in Figs. |Hl and El The MGA results with x — 1 shown in Figs. [7(a), [7(b), ED(a), 
and|^a) are very similar to the EMA results in Figs.[Ha), (He), 12(a), andEJa), respectively; 
indeed, the results for Jimd and K\ e are almost identical in the two approximations. For 
the larger value of Kxa/K\b, the enhancement factor pb&b is smaller in Fig. UJb) than in 
Fig.^c). In general, the results for this version of MGA do not show the dramatically large 
increases in Jimd seen in Fig. El for the larger ratio of Kia/Kib- In Figs. UJc), d(d) , E(b), 
and EJb) we show the same MGA calculations for pa^a, Pb^b, Jimd, and K le but with 
x = 3. The behavior does not differ greatly from that seen in the cases with x — 1, except 
for an increase in the magnitudes of Jimd and K\ e for the larger ratio of Kia/Kib- By 
contrast, the MGA results for a B host show some dramatic peaks for the larger ratio of 
Kia/Kib, as shown earlier. 



IV. DISCUSSION 



In this paper, we have calculated the intermodulation critical current Jimd i n an inhomo- 
geneous 2D superconductor characterized by a binary distribution of energy gaps. To carry 
out this calculation, we used an analogy between the effective cubic nonlinear response of an 
inhomogeneous superconductor and the effective cubic nonlinear susceptibility of a compos- 
ite dielectric medium. Using this analogy, we can apply the formalism previously developed 
to treat the nonlinear dielectric composite to the inhomogeneous superconductor. We found 
that the cubic nonlinear response of the superconductor could be expressed in terms of the 
cubic response of each "component" (i.e., energy gap), and two enhancement factors, each 
describing the field and current distribution in the related linear medium. 

In order to simplify our calculations, we have assumed that the superconducting layer 
has a binary distribution of energy gaps, and A B (with A B > A a), and we have 
considered three plausible topologies: "effective-medium" topology (A and B symmetrically 
distributed), and two "Maxwell-Garnett" topologies (A embedded in B and B embedded in 
A). We have treated all three using a simple nonlinear decoupling approximation. 

The results for Jimd are dramatically dependent on the assumed topologies. The EMA 
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and the MGA with B in A give rather similar results for ratios of Ag/A^ close to unity, 
and only modest enhancements of Jimd at any concentrations of A. However, the MGA 
with A (the component with the smaller gap) embedded in B leads to huge enhancements 
in Jimd compared to its value in either pure A or pure B, provided that the two gap values 
are sufficiently different and that x is large. 

In view of these differences, it is of interest to compare our results with the detailed 
measurements of Davis et al. [10]. These experiments do not provide results directly for 
Jimd- However, they do provide some hints about a possible connection. In particular, 
in experiments on as-grown Ni-doped samples, Ni scattering resonances were observed only 
in the regions where A < 50 meV, which were identified as superconducting regions, i.e., 
regions with superconducting phase coherence. In our results, most of the enhancement in 
Jimd comes from the enhancement factor e#, which corresponds to the component with 
a low energy gap. Therefore the regions of enhanced Jimd correspond to regions of small 
energy gap, and also regions of enhanced phase coherence according to the measurements 
of Ref. [loj. 

A striking feature of our results is the large difference between the EMA results and 
the MGA results, especially for binary composites with a large x and very different energy 
gaps. Which of these approximations is the most correct? In fact, there is not a single 
correct answer for all materials: the correct choice depends on the actual topology of the 
material of interest. In particular, we do not know, beforehand, whether the energy gaps in 
an inhomogeneous superconductor are distributed at random throughout the Cu02 planes 
or whether regions with one magnitude of energy gap are preferentially surrounded by those 
of the other energy gap. This topology determines whether we should use the EMA or the 
MGA approach. 

In the experimental gap maps the low-A regions are surrounded by the high-A 
regions in the underdoped BSCCO-2212 sample, but the high-A regions are surrounded by 
the low-A regions in the slightly overdoped as-grown BSCCO-2212 sample. Therefore, it 
appears that we can apply the MGA method to both cases, but with the roles of inclusion 
and host reversed in each case. But even this description of the distribution may be a 
simplification of the true gap distribution, which is probably continuous, not a discrete 
binary distribution. Ideally, we should consider such a continuous distribution of energy 
gaps. 
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The great sensitivity of Jimd to composite geometry, as found in the present work, is not 
surprising, in view of earlier work on transport in linear and nonlinear composite conductors 
and dielectric media. For example, it is well known that the critical exponents describing 
transport, especially nonlinear transport, in composite media are sensitive to features of the 
local geometry 20, 21, 1^. We speculate that, depending on the precise nature of this 
geometry, Jimd either diverges or goes to zero near a percolation threshold according to an 
appropriate critical exponent. 

In summary, we have presented a general formalism for calculating the intermodulation 
coefficient, and the corresponding intermodulation supercurrent density Jimd, of an inhomo- 
geneous superconductor. We have also given a simple way to calculate Jimd approximately 
in several geometries. Since such inhomogeneities are known to exist in many of the high-T c 
cuprate superconductors, this formalism is directly relevant for treating an important prop- 
erty of these materials. We find that the resulting Jimd is very sensitive to the exact spatial 
distribution of gaps within the inhomogeneous layer and thus may increase or decrease, 
depending on the topology. Our calculations show that one way to achieve a large Jimd 
is to have the component with the smaller gap and larger superfluid density embedded in 
the component with the opposite properties. This appears to be the topology seen in the 
underdoped BSCCO-2212, which thus may be well suited for a material with a large Jimd- 

Finally, we comment briefly on possible device implications of the present results. A 
useful microstrip resonator would usually have a minimum of intermodulation, or interfer- 
ence between different frequencies j^J. To achieve this, one would probably want a Jimd 
which is as large as possible, because this would lead to 1/A 2 (T, J) which has the weakest 
dependence on current. Surprisingly, we find here that Jimd can actually be increased in 
some inhomogeneous superconductors, provided that the topology is suitable. Thus the 
inhomogeneity which appears to be unavoidable in the high-T c cuprates may even be an 
advantage in constructing useful microwave devices. 
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FIGURES 




FIG. 1: Enhancement factors pa&a andpses for a 2D inhomogeneous superconductor consisting of 
a binary composite with two different energy gaps, which have the ratio K\a/K\b- The effective- 
medium approximation (EMA) is used in (a) and (c), while the Maxwell-Garnett approximation 
(MGA) is used in (b) and (d), assuming that A is surrounded by B. We use two different sets of 
ifi's: K 1A = 6.28xlO n esu/(cm-s), K 1B = 2.42x 10 n esu/(cm-s) and K 1A = 3.49x 10 n esu/(cm-s), 
K\b = 2.42 x 10 11 esu/(cm • s) as indicated in the legend. In this and the following plots, the unit 
of K\a and K\b is 10 11 esu/(cm • s). 
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FIG. 2: Intermodulation critical supercurrent density Jimd for the 2D inhomogeneous supercon- 
ductor shown in Fig. ^ The EMA method is used in (a), while the MGA method is used in (b), 
assuming that A is surrounded by B. The two sets of K\s are the same as in Fig. ^ Note 10 18 
power scale on the y-axis. 
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FIG. 3: Effective superfluid density K\ e for a 2D composite having the same values of K\a and 
K\b used in the calculations of Figs. ^ and HI (a) EMA method; (b) MGA method, taking B (the 
component with the smaller superfluid density) as the host. 
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FIG. 4: Same as Fig.^ except that x = 3. The two sets of K\s are K\a = 20.35 x 10 11 esu/(cm-s), 
Kib = 1.16 x 10 11 esu/(cm • s) and K 1A = 3.49 x 10 11 esu/(cm • s), K 1B = 1.16 x 10 11 esu/(cm • s). 
The full curve is scaled to the right axis while the dashed line to the left in (d). 
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FIG. 5: Same as Fig. ^ except that x = 3; so the two sets of K\s used are K\a = 20.35 x 
10 11 esu/(cm • s), K XB = 1.16 x 10 11 esu/(cm • s) and K 1A = 3.49 x 10 11 esu/(cm • s), K 1B = 
1.16 x 10 11 esu/(cm • s). The full curve is scaled to the right axis while the dashed line to the left 
in (b). 
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FIG. 7: Same as Fig. ^ except that the MGA method is used in all cases, assuming that B is 
surrounded by A. The two sets of K\ S are K 1A = 6.28xlO n esu/(cm-s), K 1B = 2.42x 10 n esu/(cm- 
s) and K\a = 3.49 x 10 11 esu/ (cm • s), K\b = 2.42 x 10 11 esu/ (cm • s) in (a) and (b), while they are 
K 1A = 20.35 x 10 11 esu/(cm • s), K 1B = 1.16 x 10 11 esu/(cm • s) and K lA = 3.49 x 10 11 esu/(cm • s), 
K\b = 1.16 x 10 11 esu/(cm • s) in (c) and (d). 
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FIG. 8: Intermodulation critical supercurrent density Jimd for the 2D inhomogeneous supercon- 
ductor shown in Fig. [7| The MGA method is used in (a) and (b), assuming that B is surrounded 
by A. The two sets of K\s used are the same as in Fig. Note 10 18 power scale on the y-axis. 
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